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Abstract 

A new way of magnetization switching employing both the spin-transfer torque and the torque by a magnetic field 
is proposed. The solution of the Landau-Lifshitz-Gilbert equation shows that the dynamics of the magnetization 
in the initial stage of the switching is similar to that in the precessional switching, while that in the final stage is 
rather similar to the relaxing switching. We call the present method the relaxing-precessional switching. It offers a 
faster and lower-power-consuming way of switching than the relaxing switching and a more controllable way than 
the precessional switching. 
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1. Introduction 

The high-speed and low-power-consuming mem- 
ories and data storage devices are expected in the 
current days. One of key technologies realizing 
such devices is the manipulation of the magneti- 
zation of a nanomagnet by use of the spin-transfer 
torque, which is proposed by Slonczewski [1] and 
Berger [2] . There, a spin-polarized current is intro- 
duced into an assembly of nanomagnets. The spin- 
polarized current, whose polarization direction is 
p, acts on a nanomagnet whose magnetization 
direction vector is rh. The torque caused by the 
current is proportional to the vector m x (p x rh) 
[1]. It fluctuates the magnetization [3] or even 
reverse its direction [4] . 

To date, two different methods of the magneti- 
zation switching utilizing the spin-transfer torque 
have been proposed. The first method, which 
we call the relaxing switching, is more common. 



There, the spin-polarization of the current which 
acts on the nanomagnet is almost parallel (or 
antiparallel) to the magnetization. The magne- 
tization switches its direction according to the 
directions of the spin-polarization and the flow 
of the electrons. The switching dynamics is un- 
derstood as a result of a relaxing process toward 
the equilibrium state. The main problem in the 
method for the application to memories is that it 
requires relatively high current density. The crit- 
ical dc current, Jc, is of the order of lO^A/cm^. 
A naive explanation for the need of high current 
density is the following. The magnitude of the 
torque is small when the relative angle 9 between 
the magnetization of the nanomagnet and the 
spin-polarization of the conduction electrons is 
small, since |m x (p x m)| ~ sin 6* w 6* ^ 1. 

Another problem is that the switching is slow. 
The results of the simulation indicates that the 
magnetization repeats the precessional motion un- 



Prcprint submitted to Elsevier Scienee 



6 February 2008 



til it relaxes to the equilibrium state, demanding 
a time of the order of 1 ns in the case where the 
current is 2-3 Jc [5,6]. An experimental result show 
that at least four times larger current than the 
calculated critical current is needed to switch the 
magnetization in 100 ps [7]. In summary, the relax- 
ing switching does not bring out the full potential 
ability of the spin-transfer torque. 

The second method, which are called the pre- 
cessional switching, has a possibility of overcom- 
ing such disadvantages [8,9]. There, the spin- 
polarization of the current is perpendicular to 
the magnetization of the nanomagnet. By taking 
such a configuration, it is possible to transmit the 
spin-transfer torque to the nanomagnet efficiently, 
and therefore realize the fast and low-power- 
consuming magnetization switching. However, it 
requires high-level controllability of the shape of 
the current pulse, including the pulse amplitude 
and the temporal width. Furthermore, the final 
state of the magnetization depends on its initial 
state, which indicates that the read-before-write 
is necessary. These disadvantages are due to the 
fact that the precessional switching employs the 
nonequilibrium state of the magnetization rather 
than the relaxation to the equilibrium state. 

In the present work, an alternative method, 
which resolves the problems of the previous meth- 
ods, is proposed. We call this new method the 
relaxing-precessional switching. It is different from 
the previous methods since it requires the simul- 
taneous introduction of the spin-polarized current 
and the external magnetic field. The stable mag- 
netization states under the spin-polarized current 
and the external magnetic field have been exhaus- 
tively studied in the previous work [10]. 



2. Phase diagram based on the linear 
stability 

We consider a system composed of two ferro- 
magnetic layers separated by a nonmagnetic layer. 
When the current in the direction perpendicular 
to the plane is introduced, the magnetization of 
the thinner (free) layer receives torques from the 
current, which is spin-polarized in the direction p. 



Here, p denotes the magnetization direction vector 
of the thicker (pinned) layer. In addition, an exter- 
nal magnetic field ifext is applied in the in-plane 
hard axis direction. The magnetization direction 
vector m = M/M of the free layer obeys the fol- 
lowing LLG equation: 



dm , dm 

-—— = — am X — — h 7m X ii eff 

dt dt ' 

+ 2-KM^jm X (p X rh). 



(1) 



where 



1 2/X_B gle 



2nM^ VM \e\ ' 



(2) 



Here, Hk, ol, 7, /is, y, /g, and g denote the 
anisotropy field, the damping constant, the gyro- 
magnetic ratio, the Bohr magneton, the volume 
of the free layer, the electric current, and the ef- 
ficiency factor[l], respectively. The directions of 
Jfoxt and p are chosen so that iJext = He:x.tey and 
p = e-x respectively in the relaxing-precessional 
switching. The linear stability condition for the 
equilibrium states of the magnetization [10] allows 
one to draw the phase diagram for the stable mag- 
netization states, which is shown in Fig. 1. We no- 
tice that there are bistable (B), monostable (M), 
and unstable (U) regions. In particular, we pay a 
special attention to the transition from a bistable 
region to a monostable region since it can be used 
as a switching mechanism. When the system un- 
dergoes such a transition by introducing the mag- 
netic field and the spin-polarized current, the final 
direction of the magnetization is uniquely deter- 
mined regardless of the initial direction. This gives 
the relaxing-precessional switching an advantage 
over the precessional switching. In addition, Fig.l 
indicates that the direction of the magnetization 
can be changed to desirable direction between the 
two by choosing the direction of the spin-polarized 
current. 

The critical current j^^, needed for the reversal, 
is given by the boundary line between the bistable 
region and the monostable region. The monostable 
region appears because an equilibrium state, which 
was labeled as (S) in Ref.[10], become unstable 
above the boundary line. The analytical expression 
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for this line is obtained by solving the equation (see 
Appendix for the derivation): 

i?ext = Hk ■ h{jfP ■ 2TrM/HK), 
where h{^) is a function defined by 

and X{£^) is a root of the equation : - 



(3) 



(4) 



1) = 0. The asymptotic form of around Hcxt = 
and that around i?cxt — Hk are useful. They are 
written as 

Hk 



J 



(ilext « 0) = 
1 ^^oxt 

'V2~H^ 



■ X 
27rM 

1 / Hcxt 



Hk 



O 



H, 



Hk 



(5) 




respectively. 

We notice that the critical current for the pre- 
cessional switching [9,10] = Hk/4:TtM coin- 
cides with j^^(i?cxt — 0) by setting iJoxt — 
in eq.(5). As the field i?cxt increases from zero, 
the critical current j^^(i?oxt) decreases from this 
value. Taking -ffext 0.5Hk, for example, we ob- 



tain j 



RP _ 



0.37j,f . Furthermore, by taking a = 



0.01 and Hk/4.ttM = ISOOe/l.ST = 0.0083 as 
typical values, the critical currents of the relaxing- 
precessional, relaxing, and precessional switching 



are j, 



RP 



0.003, jR(= a) = 0.01, jl 



0.008, 



respectively, which implies that the relaxing- 
precessional switching is the lowest-power con- 
suming way among the three methods. 



3. Numerical calculation of the LLG 
equation 

As was discussed in the previous work [10], the 
final state which the magnetization obtains is be- 
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Fig. 1. Phase diagram of stable states based on 
the hnear analysis. The horizontal and the vertical 
axes denote the dimensionless magnetic field J?ext/27rAf 
and the dimensionless current j, respectively. We set 
Hk/^ttM = 150Oe/0.9T = 0.017. Each of the arrows and 
the circles with dots represents the direction of a stable 
magnetization state. The former (latter) corresponds to 
the direction in (out of) the yz-plane. The labels "B", 
"M", and "U" denote bistable, monostable, and unstable 
regions, respectively. 



yond the framework of the linear stability. There- 
fore, the numerical calculation of the LLG equa- 
tion is necessary. In addition, switching speed is 
also discussed from the results. Here, it is assumed 
that the magnetic field is applied synchronously 
with the spin-polarized current. 

Fig. 2 (a) shows a time-dependent solution of the 
LLG equation on the basis of the macrospin (single 
domain) model. In this example, the parameters 
are taken as iJcxt = 0.5Hk and j = 0.003. It is 
worth comparing this solution with (b)that in the 
relaxing switching or (c)that in the precessional 
switching. 

In the initial stage of the rclaxing-precessional 
switching, the magnetization dislocates its direc- 
tion from the initial one and begins precessional 
motion in response to the introduction of the spin- 
polarized current. This dynamic behavior is similar 
to that in the precessional switching. In the typ- 
ical example of the relaxing-precessional switch- 
ing shown in Fig. 2(a), the reversal time is approx- 
imately 0.2 ns, which is almost the same as that 
in the precessional switching with a larger current 
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Fig. 2. The time evolution of mz in (a) relaxing-precessional 
switching with Hcxt = 0.5Hj^ and j = 0.003, (b) relaxing 
switching with j = 0.03, and (c) precessional switching 
with j = 0.01. We set Hk/2itM = 0.017. 



[see Fig.2(c)]. 

Once the magnetization vector gets close to ttie 
opposite direction, it relaxes to the destination. 
The dynamics in this stage is very much different 
from that in the precessional switching and rather 
similar to that in the relaxing switching. This phe- 
nomenon is advantageous for the application to 
memories since it implies that any strict control of 
the temporal width of the current pulse is unnec- 
essary. 
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Fig. 3. Phase diagram of the dynamical behavior of the 
magnetization. We set Hk/2ttM = 0.017. The definition 
of the axes is the same as that in Fig.l. The solid lines rep- 
resent the boundary lines which distinguish the final states 
of the magnetization. In "S", the magnetization stays in 
the neighborhood of the initial direction. In "P" , the mag- 
netization continues a precessional motion. In "R", the 
magnetization reverses its direction. The dashed line rep- 
resents the boundary line between the monostable region 
and the bistable one, which have been shown in Fig.l. 

4. Ranges of the current and the magnetic 
field 



The next task is to reveal the ranges of the spin- 
polarized current and the magnetic field which 
can be used as the relaxing-precessional switching. 
Here, we assume that the magnetic field is applied 
constantly in time. 

Fig. 3 is a phase diagram which illustrates the 
type of the dynamical behavior of the magnetiza- 
tion after the long-time introduction of the given 
magnitude of the spin-polarized current and that of 
the magnetic field. We have obtained this phase di- 
agram by performing numerical calculations of the 
macrospin LLG equation. In the simulation, the 
current is assumed to rise instantaneously and be 
kept a constant value since then. The whole region 
is divided into one of the three regions labeled by 
"S" , "R" , and "P" , respectively. In the case where 
the point corresponding to the given values of the 
spin-polarized current and the magnetic field lies in 
the region S, the magnetization stays in the neigh- 
borhood of the initial direction during the intro- 
duction of the spin-polarized current. In the case 
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where it lies in the region R, the magnetization 
reverses its direction and does not return to the 
initial one after the spin-polarized current is intro- 
duced. In the case where it lies in the region P, the 
magnetization precesses and approaches neither 
the initial direction nor the opposite one. There- 
fore, the region R denotes the ranges of the spin- 
polarized current and the magnetic field which can 
be used as the relaxing-precessional switching. The 
dashed line represents the boundary between the 
monostablc and the bistable states, which have 
been shown in Fig.l. 

We have two remarks. First, the monostable re- 
gion in Fig.i docs not overlap with the region S in 
Fig. 3. This implies that the monostability is a suf- 
ficient condition of the switching. Thus, the ana- 
lytical expression of the critical current obtained in 
Sec. 2 overestimates the actual switching current. 

Second, the ranges of the spin-polarized current 
and the magnetic field for the relaxing-precessional 
switching are limited to the triangle-shaped region 
shown as darkened one in Fig. 3. When the spin- 
transfer torque is large, the switching is realized 
as a precessional switching rather than a relaxing- 
precessional switching. Therefore, too large spin- 
polarized currents are not suitable. Too large mag- 
netic field is also undesirable for practical use. The 
magnetic field can be applied either by a current 
field or by a permanent magnet. In the former case, 
an additional power-consumption besides the cur- 
rent flowing through the multilayer must be taken 
into account while no additional power is consumed 
in the latter case. 

To design a practical device including the read- 
out function, we can add a magnetic reference 
layer, whose magnetization direction is collinear to 
that of the free layer, to the system. In this system, 
the effect of the reference layer to the spin-transfer 
torque is negligible since the magnitude of the 
current used in the relaxing-precessional switch- 
ing is too small to cause the relaxing switching 
as explained in the previous sections. Therefore, 
the adequate ranges of the magnetic field and the 
spin-polarized current discussed in this section are 
not affected by the addition of the reference layer. 



5. Summary 

In summary, we have proposed the relaxing- 
precessional magnetization switching, which uti- 
lizes both the spin-transfer torque and torque by a 
magnetic field. The present method can reverse a 
magnetization one order of magnitude faster with 
a current one order of magnitude lower than the 
conventional relaxing switching. In addition, it 
does not require strict control of the current pulse 
unlike the precessional switching. 



Appendix A. Derivation of eq. (3) 

In this appendix, we present the derivation of 
the analytical expression (3) of the boundary line 
between the region M and B. We concentrate on 
the case where j > and < i?ext < Hk -C 2ttM. 
First, we parameterize the vector m, which is on a 
unit sphere, in terms of the spherical coordinates 
9 e [0, tt] and (j) G [0, 2ir) as 

rux = sin cos (j), my = cos 6, niz = —sm6 sin (f). 

(A.l) 

The LLG equation (1) is rewritten as 

/A ^fm6)\^^f 1 o\ //A 

\g{0,^)) l + \-a ly [hj 
fi = cos (f){—4:TrM sin sin (p — Hk sin sin (j) 

+ 2nMj cos 0) 
/a = {Hk sin^ 4> - AttM cos^ (p) sin cos 

- iJext sin - 2iTMj sin ^ (A.2) 

in the spherical coordinate. The equilibrium 

states are given by the solutions of the equations: 
f{0, 0) = g{0, (j)) = 0. The solution which becomes 
unstable for large enough j > and ifext is labeled 
as (S) in Ref.[10], which continuously changes into 
m = —By as ffext < —Hk and j 0. The 

solution (S) is obtained by 

Hk sin cos ^— ilext sin 0—2TrM\j\ = 0, sin = sgn(j). 

(A.3) 

In the following, we consider the existence of this 
solution and the stability of the state correspond- 
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Fig. A.l. L-i : y = ±\/l - x-'^ and L2 : 
y = (2iTMi/Hx)x + Hgy^t / H X ■ When the two lines have 
common points, the equilibrium states corresponding to 
eq.(A.3) 
exist. 

ing to this solution. The equation (A. 3) is equiva- 
lent to 



±yjl-x-^ = {2nMj/HK)x + H,^t/HK, (A.4) 

if we set a; = 1/ sin 6. Thus, the existence of 
the solution is guaranteed when the two lines 
in the xy-planc, Li : y = ±\/l — x~'^ and L2 : 
y = {2nMj/HK)x + Hc^i/Hk, share at least one 
point. Reminding that the gradient of the line L2 
is positive, we see from Fig. A.l that the two lines 
intersect when the y-intcrccpt of the line L2 is 
smaller than some critical value h(2'!T M j / H k) , 
which is a function of the gradient of L2- In other 
words, the existence condition is written as 



ext 



< h 



K 



2nMj 
Hk 



(A.5) 



In the inequality (A.5), the equality holds when L\ 
and L2 are tangent to each other. In that case, the 
tangent point {xo,yo) = (1/ sin ^0, cos ^0) satisfies 



yo 



ax 

d f27TMj , f2TTMj\ 



1 



2nMj 



H 



-Xq 



K 



271 Mj 



H 



K 



(A.6) 



(A.7) 



X=Xq 



The condition (A.7) yields 

which is rewritten as 
2'KMj 



K 



1)=0, 



sin^ ^0 
cos^o 



(A.6 



(A.9) 



By introducing the function X{^), which is defined 
as a root of the equation : + ^^(X — 1) = 0, we 
have x, 



- X{2TTMj/HK) and 
27rAfj 



Hk 

1 - 2X 



f 2TTMj\ 



)) 



When 7?cxt < Hk ■h{2TrMj / Hk), Li and L2 inter- 
sect each other at two points. Therefore, eq. (A. 3) 
have two solutions 0< and 6>, which satisfy < 
9< <0a <9y < tt/2. 

In general, the equilibrium state is stable if all 
eigenvalues of the matrix: 



L : 



/ df/de df/dc^ 
\dg/de dgld4>j 



(A.11) 



/=s=o 

have negative real parts, in other words, tri < 
and dot L > 0. In the present system, the trace of 

L is negative because 



1 

-trL = - 2 
a 



H 



K 



2ttM 



(3cos2 6'- 2) - 2 



Hp- 



2itM 



cos 9 



< 



2+^+2 



H, 



< 0. 



(A.12) 



2wM 2itM 
Here, we have used the relation (A. 3). 

In contrast, the sign of the determinant of L 
cannot be determined automatically. Since det L is 
calculated as 

Hk 



det L ={1 + 0^) [2 + 
Hk 



2ttM 



2ttM 
.cos 6* 



■ sin^ 



HsKt 

2ttM 



cos X 



sin^ 



(A.13) 



we obtain 



sgn(deti) =sgn( 
=sgn( 



sin'^ 9 2TrMj 



cos 6 
sin^e 



Hk 
sin^ 60 
cos 6 cos 60 



) 

)=±1 



(A.14) 
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Thus, the equihbrium state corresponding to the 
solution 9y is always stable if it exists. In contrast, 
that corresponding to 0< is always unstable. 

In summary, we have proved that there exist 
two equilibrium states under the condition ffoxt ^ 
Hk ■ h{2TTMj / Hk)- One of the equilibrium states 
is always stable if it exists. Namely, the above con- 
dition is not only the existence condition of a equi- 
librium state but also its stability condition. 
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